München, 01. Juli 2002 ISSN 0341-8391 


SPIXIANA 183-191 


Tadpole morphology of two spadefoot toads 
(Pelobates fuscus and P. syriacus) 


(Amphibia, Anura, Pelobatidae) 
Vesna Sidorovska, Katarina Ljubisavljevic, Georg Dzukic & Milos L. Kalezic 


Sidorovska, V., K. Ljubisavljevic, G. Dzukic & M. L. Kalezic (2002): Tadpole 
morphology of two spadefoot toads (Pelobates fuscus and P. syriacus) (Amphib- 
ia, Anura, Pelobatidae). — Spixiana 25/2: 183-191 


Intraspecific variation in body mass and morphometric and qualitative char- 
acters of different larval stages were studied in a number of population samples 
of Pelobates fuscus and P. syriacus from Serbia and Macedonia. Tadpoles were 
collected from populations within the areas of allopatry and sympatry, in both 
allotopic and syntopic locations. It was discovered that, on the average, larval 
body size increased linearly to stage 41, and then suddenly decreased in the last 
stages. Changes in size and shape relationships, studied by Principal Compo- 
nent Analysis, revealed that interspecific differences appeared to be much 
larger between younger tadpoles, with the main difference being in head size 
measurements. Convergence in morphometric aspect progressed towards the 
time of metamorphosis. Canonical Discriminant Analysis revealed slight spe- 
cies differences on the basis of tadpole morphometric characteristics, including 
body mass. This was found for both pooled samples of all stages analysed, as 
well as for particular larval stages. However, on the whole, the effect of locality 
was apparently very influential. Studied qualitative traits failed to discern any 
geographic trends or provide a ready means to segregate species among the 
different localities. Consequently, it can be inferred that tadpole morphological 
characters can not be used to distinguish between the two species except in 
early growth stages. 
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Introduction 


The spadefoot toads of the genus Pelobates in- 
clude four species, among them the common 
spadefoot (P. fuscus) and the eastern spadefoot 
(P. syriacus). These two toads are clearly differ- 
entiated species. They belong to two different 
phylogenetic lineages within the genus Pelo- 
bates (Maglia 1998). With respect to their distri- 
bution, P. fuscus has an extensive range which 
covers most of the plains and hilly steppic re- 
gions of Central, Eastern and Southeastern Eu- 
rope (Nollert 1997). The eastern spadefoot 
(P. syriacus) has a range extending from the 
southern part of the Panonian Plain (Banat) 
through Southeastern Romania, Bulgaria, and 
Greece to Southwestern Asia (Sofianidou 1997). 
P. fuscus and P. syriacus have narrow zones of 
sympatry. Today these two species occur to- 
gether along the lower course of the Danube, as 
well as in the vicinity of the Bosphorus to the 
south (Dzukic & Pasuljevic 1983, Eiselt 1988). 

Adult specimens of P. fuscus and P. syriacus 
can be easily distinguished due to pronounced 
differences in coloration patterns, as well as by 
the presence (P. fuscus) or absence (P. syriacus) 
of a well-marked dome on the top of the head. 
Also, adult individuals of P. syriacus are on the 
average, much larger than those of P. fuscus, 
including within the area of their sympatry 
(Rot-Nikcevic et al. 2001). These two species 
also differ in the direction of the magnitude of 
sexual size dimorphism. P. fuscus follows the 
general trend in amphibian species, having 
much larger females than males, while in 
P. syriacus sexes have approximately the same 
body size, or the males can be larger than fe- 
males in some cases (Rot-Nikcevic et al. 2001). 
On the other hand, studies of life-history traits 
showed that interspecific differences in longev- 
ity and time of attainment of sexual maturity 
appeared to be small and without a consistent 
pattern of variation. 

Tadpoles of all European anuran congener- 
ic species can be distinguished by their external 
morphological traits, except for green frogs 
(R. esculenta complex) and spadefoot toads. The 
search for possible diagnostic morphological 
traits for their tadpoles is of importance in the 
area of species sympatry, especially in the sites 
of syntopy (i.e. where adults of both species 
breed concurrently). 
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To the best of our knowledge, an in-depth 
morphological analysis of Pelobates tadpoles, 
especially intra- and interpopulation variation 
patterns, has never been done. Current knowl- 
edge comes primarily from a few publications 
describing only a limited number of tadpole 
traits, for example body shape and coloration. 
In regard to interspecific comparisons, in sev- 
eral taxonomic keys it was claimed that tad- 
poles of P. fuscus and P. syriacus can not be 
distinguished (e.g., Engelmann et al. 1986), 
while Fuhn (1960) and Kuzmin (1999) pointed 
out some differences, mainly in coloration and 
body shape. 

The purpose of this paper is twofold. First, 
an in-depth analysis of both morphometric and 
qualitative external morphological characteris- 
tics of different larval stages in a number of 
population samples of P. fuscus and P. syriacus 
was performed. Second, in analysing tadpoles 
which came from the areas of allopatry and 
sympatry, both allotopic and syntopic locations, 
this study attempts to evaluate the influence of 
location, stage of tadpole development and 
patterns of data variation for both species. 


Material and methods 


Population samples. Tadpoles were collected 
from allopatric areas: three localities for P. fus- 
cus (Hrastovaca, Crepaja and Samos) and two 
localities for P. syriacus (Prdejci and Djavato), 
as well as from the zone of sympatry (Deliblato 
Sand). Within the range of species coexistence, 
the allotopic breeding sites of P. fuscus, but not 
of P. syriacus, and more numerous syntopic 
breeding sites were recognized. In this analy- 
sis, one allotopic population sample (Zamfir 
bara) and two nearby population samples with 
tadpoles of both species (Zubanov and Djabin 
salas), have been included (see Appendix for 
sample size and more details about localities). 
Samples were collected in June and July of 1999 
and 2000, and preserved in 70 % ethanol. Ana- 
lysed larvae were in various stages of develop- 
ment: from 34 to 45 (according to Gosner 1960). 
In order to obtain an appreciable number of 
larvae for our analyses, we pooled initial stages 
into reduced stage orders as follows: stage I 
(stages 34, 35 and 36); stage II (stages 37 and 
38); stage III (stages 39 and 40); stage IV (stage 


Fig. 1. Body size characters measured. For character abbreviations see text. 


41); and stage V (stages 42, 43, 44 and 45). 
Larvae from Djavato were kept in the laborato- 
ry for a period of three months (March-June, 
2000), from the blastula to tadpole stage 44. 
Larvae were housed in an aquarium (150 1), 
with various plants from their natural habitat. 
The aquarium was exposed to a natural pho- 
toperiod and daily changes in room tempera- 
ture (15-25 °C). A constant water level was main- 
tained. Water was continually recycled and fil- 
tered. Thus, the larvae were maintained in 


standard conditions with ad libitum access to 
food (boiled spinach and commercial aquari- 
um fish food). At the beginning of June, when 
larvae had attained a size between stages 36 to 
44, they were preserved by immersion in 70 % 
ethanol. 

The specimens are deposited in Georg Dzu- 
kic’s Batrachological collection (Institute for 
Biological Research, Belgrade), and in the Ver- 
tebrate collection of the Macedonian Museum 
of Natural History, Skopje, Macedonia. 
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Characters studied. We measured the follow- 
ing characters with a digital caliper to the near- 
est 0.01 mm (Fig. 1): Drn: snout to nostril dis- 
tance (distance between the tip of the snout to 
the nasal pore), Dpno: posterior nostrils to an- 
terior eyelid distance (distance between the 
anterior eyelid commissures to the posterior 
edge of the nasal pore), Drs: snout to spiracle 
distance (distance between the tip of the snout 
to the spiracle pore), Dspp: spiracle to hind- 
limbs distance (distance between the spiracle 
pore to the anterior margin of the hindlimb at 
insertion), Altmcd: tail musculature height (at 
insertion of hindlimbs), Lep: body length (from 
the snout to the anterior margin of the hind- 
limb at insertion), Led: tail length (from the 
anterior margin of the hindlimb at insertion to 
the tail tip), L: total length (from the snout to 
the tail tip), Altcd: tail height (maximum tail 
height), Dro: snout to eye distance (distance 
between the tip of the snout to the anterior 
eyelid commissures), Dano: anterior nostrils to 
anterior eyelid distance (distance between the 
anterior eyelid commissures to the anterior edge 
of the nasal pore), Dra: snout to anal pore dis- 
tance (distance between the tip of the snout to 
the posterior edge of the anal pore), Dsa: spir- 
acle to anal pore distance (distance between the 
spiracle pore to the posterior edge of the anal 
pore), Ltro: width of mouth region, Spi: dis- 
tance between nostrils, Spp: minimum interor- 
bital distance, Spo: interocular distance (dis- 
tance between the eyes), Lo: eyeball length 
(from the anterior to the posterior eyelid com- 
missure), Ltcn: head width at position of nos- 
trils, Ltco: head width at position of posterior 
eyelid commissure, Ltc: maximum head width, 
F: length of femur (from insertion to the distal 
end of the femur), T: length of tibia (from the 
proximal end of the tibia to the metatarsal wrist), 
N: foot length from the heel to the end of the 
longest toe, P: foot length from the metatarsal 
wrist to the end of the longest toe, Ltf: width of 
femur (at the widest part), Ltt: width of tibia 
(at the widest part), Cint: length of inner met- 
atarsal tubercle. Body weight (M) was meas- 
ured to the nearest 0.001 g. 

Qualitative characters were as follows: 
I: Colour of body: a: golden yellow, b: olive 
green, c: grey, d: black; II: Presence of spots: 
a: present, b: not present; III: Shape of margin- 
al papillae of mouth region: a: conical, b: py- 
ramidal, c: low conical; IV: Constitution of tad- 


186 


pole body: a: stout, b: slender; V: Shape of tail 
tip: a: rounded, b: pointed. 


Statistical procedures. In order to minimise 
deviations from normality caused by allomet- 
ric relationships, all body measurements and 
body mass were log-transformed. Student's 
t-test was used to test differences in mean val- 
ues between localities for every stage, consid- 
ering P<0.05 as a level for significance. A one- 
way analysis of variance (ANOVA) was per- 
formed using population samples or stages as 
factors. In order to determine the degree of 
difference between species and among stages, 
a multivariate analysis of variance (MANOVA) 
was used. Principal Component Analysis (PCA), 
was used to examine patterns of morphometric 
variation between species through ontogenic 
changes in tadpoles’ size and shape. The anal- 
ysis was performed on the covariance matrix of 
log-transformed variables. Separate PC analy- 
ses were computed for five larval stages of two 
samples (Hrastovaca: P. fuscus and Prdejci: 
P. syriacus). Canonical Discriminant Analysis 
(CDA), of log-transformed data was used to 
distinguish between intraspecific and interspe- 
cific variation for morphometric measurements. 
To reduce the influence of body size, CDA was 
also performed on residuals from a regression 
analysis of each trait on snout to anal pore 
distance (Dra). Percentages of states were cal- 
culated for each qualitative trait. Programs de- 
signed and implemented by the STATISTICA 
(StatSoft, Inc. 1997) package were used. 


Results 


Morphometric characters. The largest meas- 
ured P. fuscus tadpole came from the Samos 
sample, with a body length of 110.8 mm (stage 
HI), while the minimum total length for a 
P. fuscus tadpole was recorded in the Zamfir 
bara sample (53.5 mm, stage I). Within P. syria- 
cus, the maximum and the minimum total 
lengths were recorded in the Prdejci sample 
(104.0 mm, stage IV, and 40.8 mm, stage V). 
Univariate analysis of variance demonstrated 
high significance among sample differences for 
all characters (ANOVA, P<0.001). Moreover, 
there were significant differences among sam- 
ples within species. Consequently, tadpoles 
from a pair of P. syriacus samples (Prdejci and 


Djavato) significantly differed in 17 out of 29 
characters. Among P. fuscus samples, a com- 
parison of tadpoles from Hrastovaca vs. Samos 
and Crepaja had 21 out of 29 characters being 
significantly different. On the other hand, 
P. fuscus tadpoles from two geographically close 
samples (Samos and Crepaja) differed signifi- 
cantly in only 6 of 29 traits. The greatest differ- 
ences among population samples appeared to 
be in comparisons of the first and the second 
larval stages (significant differences in 27 of 29 
characters, and in 26 of 29 characters, respec- 
tively). However, differences were small or 
without significance in the third and the fourth 
stages with regard to hindlimb measurements 
(significant differences in 1 of 7 hindlimb char- 
acters, and in 2 of 7 hindlimb characters, re- 
spectively). In the fifth stage, differences be- 
came the smallest for most morphometric char- 
acteristics (14 non-significant traits), except for 
hindlimb measurements. 

When the larval stages were considered sep- 
arately, there was a tendency for increased tad- 
pole body length (L) up to the fourth stage 
(stage 41 according to Gosner 1960), and a sub- 
sequent rather sharp decrease in the last stage 
(V —42-45) due to drastic metamorphic changes 


(Table 1). On the average, the largest tadpoles 
at the fourth stage were P. fuscus from the Zam- 
fir bara sample (100.1+5.7 mm), and the small- 
est, P. syriacus tadpoles raised in the laboratory, 
Djavato sample (72.6+1.6 mm). Taking into con- 
sideration all studied morphometric characters, 
paired testing of differences among samples 
revealed a complete departure of the Djavato 
sample (P. syriacus) at every stage. When the 
tadpoles of the first larval stage were com- 
pared, the greatest differences in pairwise sam- 
ple comparisons appeared between Djabin salas 
(syntopic site) and Djavato (significant differ- 
ences for 16 of 22 character comparisons). Non- 
significant differences were found for Zamfir 
bara (P. fuscus) vs. Zubanov salas (syntopic 
site); Zamfir bara vs. Prdejci (P. syriacus) and 
Zubanov salas vs. Djabin salas (syntopic sites) 
pairs. At the second stage the greatest differ- 
ences existed when the Djavato sample (P. syria- 
cus) was compared with the two P. fuscus sam- 
ples (Hrastovaca and Samos, 19 out of 29), while 
the smallest differences appeared between two 
nearby P. fuscus samples (Samos and Crepaja, 3 
out of 29). The greatest differences at the third 
stage were between samples from Hrastovaca 
and Crepaja (P. fuscus, 18 out of 29), while the 


Tab. 1. Sample means and standard errors (SE) for total length: L (in mm) in different tadpole stages of 


population samples. N = sample size. 


I Il TI IV y 
Locality N N N N N 
Mean +SE Mean+SE Mean+SE Mean+SE Mean + SE 
octavae fuscus) 71 swab 75.98 21.80 PPAT 85.26 Te T 7.87 
BERLE a Rca o 1.90 ee spice aean E 
Crepaja (P. fuscus) a nee ne are TA 1.95 3 
aadik pereda LE Hane L : 1001558 = 
Zubanov salas (P. fuscus/P. syriacus) 6 MR 43 p P ee 6.84 7 
e eleo e e ee, e h sae Pa E eR ae 6.02 
E ea) ee T 7.92 TA 8.42 eae Tes ep 
Djavato (P. syriacus) ae : 1.59 72.63: 3.86 a 72615163 68.875170 
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samples from Samos and Zubanov salas lacked 
any significantly different character compari- 
sons. The fourth stage revealed the greatest 
distinction between tadpoles from two P. syria- 
cus samples (Prdejci vs. Djavato, 21 out of 29), 
while only one out of 29 character comparisons 
appeared to be statistically significant when 
Zamfir bara and Zubanov salas were compared. 
The testing in the fifth stage could include only 
four samples. The greatest differences were 
between two P. syriacus samples (15 out of 25 
character comparisons), and the smallest ones 
were between the Djabin salas and Prdejci sam- 
ples (1 out of 25). 

Multivariate analysis of variance (MANO- 
VA) with species and stage as factors revealed 
significant differences in snout to nostril dis- 
tance (Drn, P<0.05), body length (Lep, P<0.01), 
snout to anal pore distance (Dra, P< 0.001), head 
width at position of nostrils (Ltcn, P<0.01), 
maximum head width (Ltc, P<0.05), and hind- 
limb measurements (F, P<0.05; T, P<0.01; N, 
P<0,05; P, P<0.05; te 0052 Ltt, P<0.01). 

PCA showed significant differences between 
P. fuscus and P. syriacus tadpoles during larval 
stages. There was obvious segregation on the 
first principal component (PC1) in the first stage. 
The body weight (M), head width at position of 
nostrils (Ltcn), maximum head width (Ltc), 
snout to nostril distance (Drn) and spiracle to 
hindlimbs distance (Dspp) were the highest- 
loading characters on the first component. Tad- 
poles of P. syriacus had higher relative values 
for these traits. Segregation in the second stage 
is less obvious, although there was a certain 
distinction along the PC2 on the basis of a great- 
er influence from head width (Ltcn, Ltc and 
Ltco) in P. syriacus tadpoles. A tendency for 
differentiation on the basis of these same char- 
acters continued during the third and fourth 
stages where Ltcn, Ltco and Ltc were the high- 
est-loading characters on the PC3. There were 
no apparent differences between species in the 
fifth stage (figures and additional statistical data 
are available from corresponding author upon 
request). 

The Canonical Discriminant Analysis (CDA), 
of morphometric characters and body weight, 
performed on pooled population samples with- 
out regard to larval stage, revealed 82.5 % of 
the total variation expressed in the first three 
canonical axes. The mean percentage of cor- 
rectly classified cases in each sample was very 
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high (89.6 %). The smallest correctness value 
(66.7 %) was for Zubanov salas (syntopic site), 


| 
| 


and the greatest (100 %) from the samples Sa- — 


mos, Crepaja (P. fuscus) and Djavato (P. syria- 
cus). The first three canonical variates appeared 
bipolar (the table of pooled within-groups cor- 
relations between variables and canonical axes 
is available from corresponding author upon 
request). Hindlimb measurements were the re- 
verse of the other traits on the second canonical 
variate. Most traits had the highest and most 
significant correlation on the second variate. 
This suggests that discrimination among sam- 
ples based on this canonical variate may be 
caused by overall size. Differences between 
species were the most discernible in the region 
of the first and the second canonical axes (Fig. 2). 
Samples of P. syriacus and P. fuscus were situat- 
ed on opposite sides of the second axis. Tad- 
poles of P. syriacus from Macedonian samples 
(Prdejci and Djavato) had in general, higher 
values for hindlimb measurements in compar- 
ison to the others. On the other hand, they had 
the lowest values for some other traits (Drn, 
Dpno, Altcd, Dro, Dano, Spi and Spp). Larvae 
of P. fuscus from the Samos and Crepaja sam- 
ples, situated in the upper right quadrant, were 
differentiated from the remaining ones mainly 
on the basis of the longest and tallest tail (Led, 
Altcd) and the longest body (L), but they had 
the lowest mean values for all characters con- 
cerning hindlimbs. Of special note, one sample 
of presumably mixed P. fuscus and P. syriacus 
tadpoles (from Zubanov salas) appeared well 
in the morphospace of P. fuscus, while some 
individuals of the second sample from a syn- 
topic site (Djabin salas) were close to P. syriacus 
tadpoles from Prdejci. CDA performed on re- 
siduals, when the size effect was eliminated, 
did not reveal greater distinction between spe- 
cies. 

A separate CDA was done for the first and 
the fourth larval stages, since these two stages 
were present in all samples. In the first stage, 
the CDA revealed 79.1 % of the total variation 
expressed in the first three canonical axes. There 
were no misclassifications among samples (cor- 
rectness of 100 %). The apparent distinction of 
the Djavato sample (P. syriacus) was due to the 
low mean values for the characters M, Drn, 
Led, Dano, Dro, which made the greatest con- 
tribution to variance on the first canonical axis. 
The Hrastovaca sample was differentiated from 
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Fig. 2. Plots of canonical variable scores of the first two canonical axes (CAN) of the pooled samples of five 


tadpole stages. 


the remaining samples on the basis of the low 
mean values for the characters Dra, Dsa and 
Drs. In the fourth stage, the CDA revealed 
92.8 % of the total variation expressed in the 
first three canonical axes, with correct classifi- 
cations of 100 %. The CDA showed a high lo- 
cality effect but without any distinctive group- 
ing of analysed population samples (figures 
and additional statistical data are available from 
corresponding author upon request). 


Qualitative characters. Black body colour (Id) 
was predominant only in Crepaja sample (P. fus- 
cus), while grey (Ic) colour prevailed in other 
samples (the table of qualitative traits states 
percentages is available from corresponding 
author upon request). Tadpoles from most sam- 
ples had unspotted bodies (IIb) and a mouth 
region with low conical papillae (IIIc). Samples 
of P. fuscus from Crepaja and in mixed samples 
from a syntopic locality (Djabin salas), as well 
as the P. syriacus sample from Djavato, more 
frequently had spotted tadpoles (Ila) than oth- 
ers. Tadpoles of two P. fuscus samples (Hrasto- 
vaca and Crepaja) and one of P. syriacus (Djava- 
to) appeared to have their mouth region pre- 
dominantly, with conical papillae (Ma). In all 
samples, individuals with pyramidal papillae 
(IIIb) appeared with low frequency. Tadpoles 


with a slender body (IVb) and pointed tail (Vb), 
prevailed in most samples. Samples from one 
of the syntopic locations (Zubanov salas), and 
of P. fuscus (Zamfir bara) and P. syriacus (Prdej- 
ci) more frequently had tadpoles with stout 
bodies, while samples from the syntopic sites 
(Zubanov salas and Djabin salas) had exclu- 
sively or predominantly tadpoles with a round- 
ed tail. 

In a separate study, it was determined that 
the number and arrangement of keratodont 
rows did not differ between tadpoles of P. fus- 
cus and P. syriacus at any larval stage studied 
(unpubl.). 


Discussion 


Tadpole body size. It is commonly understood 
in batrachology that spadefoot toads have the 
largest tadpoles in comparison to other anurans. 
However, except for a few scattered estimates 
in the literature, body size data obtained for an 
appreciable number of larval individuals from 
the same breeding site are almost non-existant. 
At hatching, tadpoles of P. fuscus and P. syria- 
cus measure 4-6 mm (Kuzmin 1999). A total 
length of about 120 mm for P. fuscus individu- 
als prior to metamorphosis in the same year of 
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hatching has been reported, while the multiple 
overwintering larvae (three years old) can reach 
200 mm in total length (Grillitsch et al. 1983). 
To the best of our knowledge, the maximum 
length thus far recorded for a P. syriacus tad- 
pole is 165 mm (Kuzmin 1999). Body size data 
for Pelobates tadpoles recorded from the neigh- 
bourhood of FR Yugoslavia are also rare and 
are based on measurements of only a few indi- 
viduals. Thus, for P. fuscus tadpoles, Grillitsch 
et al. (1983) reported a total length from 80 to 
120 mm (Austria), up to 100 mm for tadpoles in 
Hungary (Dely 1967), and 95.0 mm for tad- 
poles in Romania (Fuhn 1960). P. syriacus tad- 
poles from Greece range in total size from 90 to 
115 mm (mean value 101.5 mm, Sofianidou 
1977). From Turkey a size of 137 mm has been 
recorded (Zaloglu 1964). Apparently, the body 
size for spadefoot toad tadpoles observed in 
this study is well within the range of P. fuscus 
and P. syriacus so far recorded. 


Body size and shape changes during onto- 
geny. On the basis of morphology, dissimilar- 
ities in body size are apparent between P. fus- 
cus and P. syriacus adults. Interestingly enough, 
spadefoot toads of these two species start ter- 
restrial life at approximately the same body 
size, but diverge soon during the first year of 
their juvenile phase (Rot-Nikcevic et al. 2001). 
Specifically, juveniles of P. syriacus appeared to 
have a significantly higher growth rate than 
that of P. fuscus, at least in the zone of sympa- 
try. Also, the larger adult size for some P. syri- 
acus individuals is the result of continued 
growth after the attainment of sexual maturity 
(the age at maturity ranges from two to four 
years; Rot-Nikcevic et al. 2001), which is not the 
case for P. fuscus toads. 

During the larval aquatic phase, differences 
among analysed samples of P. fuscus and P. sy- 
riacus followed the opposite trend. Changes in 
size and shape relationships, revealed that in- 
terspecific differences appeared to be much larg- 
er between younger tadpoles, attributable pre- 
dominantly to head size measurements, while 
the maximum convergence in morphometric 
respect progressed towards the time of meta- 
morphosis. 

Slight species differences on the basis of the 
tadpoles’ morphometric characteristics, includ- 
ing body mass, were found for pooled samples 
of all stages analysed, as well as for particular 
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larval stages. However, in the overall pattern of — 


| 
| 
; 
| 


sample differentiation, the effect of locality was — 


apparently very influential. Thus, two P. syria- 
cus samples in which individuals were grown 
under quite differing conditions (tadpoles from 
Prdejci were collected from nature while those 
of Djavato were raised in an aquarium) dif- 
fered substantially in spite of the close geo- 
graphical proximity of their habitats, which does 
not exclude a possible gene-flow between pop- 
ulations. Tadpoles of P. fuscus also experience a 
great site influence in the shaping of morpho- 
metrics, to the extent that no misclassifications 
among samples were observed in some cases. 
Individuals of one sample from a syntopic site, 
with tadpoles of both species, were placed main- 
ly within the P. fuscus morphospace, while some 
individuals of another syntopic sample ap- 
peared to be closer to one P. syriacus sample. 
In summary, from these results it can be 
inferred that multivariate morphometrics of tad- 
poles, as well as the study of their integumental 
qualitative traits (which failed to discern any 
geographic trends or segregation of species 
among the different localities), can not discrim- 
inate between spadefoot toad species, especial- 
ly those from the zone of sympatry. Thus it can 
not be used in the process of species identifica- 
tion, except in very early stages. The use of 
molecular techniques, such as genetic markers 
of allozyme loci, would need to be used to 
assess the specific status of premetamorphosed 
individuals of P. fuscus and P. syriacus. 
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Appendix 


Localities of population samples, their altitudes, 

UTM code (10 x 10 km’), and the sample size (num- 

bers of tadpoles in larval stages [+II+III+IV+V). 

1. Hrastovaca (Vojvodina, Serbia, 110 m above sea 
level, CS 91, 8+11+7+11+4) 

2. Samos (Vojvodina, Serbia, 120 m, DR 80, 4+6+3+ 
2+0) 

3. Crepaja (Vojvodina, Serbia, 95 m, DQ 78, 1+11+ 
12+3+0) 

4. Zamfir bara — Deliblato Sand (Vojvodina, Ser- 
bia, 150 m, EQ 16, 5+6+0+4+0) 

5. Zubanov salas — Deliblato Sand (Vojvodina, Ser- 
bia, 75 m, EQ 26, 4+0+2+3+0) 

6. Djabin salas — Deliblato Sand (Vojvodina, Ser- 
bia, 75 m, EQ 26, 13+34+11+5+3) 

7. Prdejci (FYR Macedonia, 66 m, FL 26, 5+3+4+ 
7+16) 

8. Djavato (FYR Macedonia, 66 m, FL 26, 10+4+1+ 
9+11) 
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